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FOUEVVOIU) 


A {>i'iu‘ral  throt'-diiiuMisional  comiJutalional  model  for  nonlinear  com- 
posite structures  and  materials  was  devidoped  hy  I’rototype  Development  Associate 
Inc.  (PDA)  under  Contract  Number  N0()014-7(i-C-^01(il  tor  the  Oltice  ol  Naval 
Uesearch.  Dr.  K.  D.  Stanton  servad  as  Program  Manager,  Mr.  I..  M.  Crain  was 
responsible  tor  scientitie  |)rogramming  services,  Dr.  D.  Mulville  served  as 
technical  monitor  tor  the  Navy  and  Dr.  .1.  Bueh  ot  PDA  contributed  a wealth  ot 
background  information  on  material  bc-havior.  Their  support  and  that  of  Dr.  N. 
Perrone,  who  originated  the  project,  are  gratefully  aeknowledged. 


A HS' IMPACT 

A eonipiiUitional  model  I'or  the  analysis  ol  sliuetural  and  microslruttural 
beluivioi’  in  general  solids  of  composite  material  is  presented.  Emphasis  is  placed 
on  representing  the  anomalous  material  behavior  of  composites  and  on  the  construc- 
tion of  compuUUional  motlels  with  variable  properties.  Alternative  material  models 
using  continuum  and  statistical  mechanics  were  reviewed  and  a modular  code  de- 
signed for  compatibility  with  several  different  models.  The  constituent  materials 
are  characterized  in  terms  of  those  state  variables  that  correlate  a materials 
response  such  as  effective  stress  or  strain  energy.  A parametric  cubic  repre- 
sentation is  used  for  all  state  variables,  the  solid  geometry  and  all  physical  pro- 
perties. The  associated  finite  element  extends  isoparametric  mcxleling  to  allow 
properties,  linear  or  nonlinear,  to  vary  over  the  volume  of  an  element  as  in  rosette 
material  construction.  Applications  of  tlie  model  Uj  a carbon-carbon  unit  cell,  to 
strain  singularities  and  to  the  inelastic  response  of  a graphite  bar  illustrate  its 
utility.  Good  agreement  with  triaxial  test  data  for  inelastic  strains  under  hydro- 
static pressure  is  obtained.  Numerical  results  are  computed  using  PATGHKS-Ill 
and  the  conjugate  gradient  algorithm  without  the  assembly  of  large  matrices.  This 
approach  is  tailored  for  vector  processors  and  can  reduce  the  high  cost  of  non- 
linear thrce-iii mensional  analyses. 
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inrri’ioN 


'I'lic  usi'  111'  composite  matci'ials  lor  primary  sti'ucturc  has  led  to  diMnaiuls 

lor  inlormation  on  their  thri’e-dimensional  response  simply  because'  uniaxial  and 

biaxial  loading  often  lead  to  impoi-tant  triaxial  el'leets.  'lliis  is  true'  in  the-  elastic- 

range  uherc-  edge-  c-flc-ets  can  cause-  de-lamination  anel  it  is  particularly  true  in  Ihe 

inc-lastie-  I’ange  \vhc-i-e  load  rc-disi  I'ibut ion  can  be  strongly'  dependent  on  the  mic-ro- 

slruetu)-c-  of  the-  matei'ial.  t hie  of  the  moi'e-  re-aelable  essays  on  this  subject  is 

Di'uckc-r's^  acceiunt  in  words  and  pictui'e-s  ol  how  elille-rent  microstructurc-s  ellc-t-t 

inc-lastic'  macrosc-opie  re-sponsc-.  Many  ol  tlu-sc-  inveilve  inelastic  mec-hanisms  ejuite- 

dilferent  trom  the  e-rystallographie  slip  characteristic  of  metals.  Some  graphites 

2 

e'xhibit  biaxial  sotteiiing  den  in  |>art  to  mic-rocrac-king  anel  se-vc-ral  cor ri-lativc- 
niMlcIs  * (i.e-.,  plu-nomenologie-ali  arc  acailable  that  tit  this  data  in  the  tension- 
te-nsion  i|uadrant.  c^uib  the-  oppeisite-  be-hac  ior  e-an  occur  in  the  compivssion- 
eompression  iiuadrant  toi-  another  particulate  ceimposile,  cemcrete,  where  the 
c-onfine-menl  of  mie-roc-racks  results  in  biitxial  stille-ning  and  again  ceirrelative 
meidels  ’ are  available-  that  tit  this  data.  Anomalous  inelastic-  behavior  also  eiccurs 
in  tibc-r-reintorci-d  c-ompositi-s  as  elescrib(-d  in  the  survey  paper  by  Francis  and 
Bert*’  whei  review  the  history  ot  the  tibc-r-matrix  load  transter  controversy.  This 
issue-  is  ceimplie-alc'd  by  mie-rostructural  de-tects,  partie-ularly  in  compositc-s 
subjc-cted  to  seve-rc-  processing  conelitions,  bv  local  tiber  buckling  and  by  residual 
preice-ssing  strains. 

Genuinely  predictive  inelastic  models  tor  composite  materials  are  not 
available  anel  most  correlative  models  arc  an  extrapolation  of  uniaxial  data  te>  multi- 
axial  stress  states.  The  situation  in  metals  is  similar'  but  not  to  the  same  degree. 
Lin  and  llavner'  , for  example,  have  prexluceel  m<Klels  that  can  predict  multiaxial 
stress-strain  respe>nse  for  polycrystalline  aggregates  and  microstructural  models 
at  the  atomic  level^**  (lattice  statics)  are  alsei  used  in  the  material  sciences.  The 
microstructural  detail  required  in  a predictive  material  mixlcl  for  inelastic  response 
clearly  can  be  prohibitive  for  structural  mechanics.  In  the  case  of  conii)osite 
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matii-i:i!s  r\ni  thr  a.-'Siiiiipl  inn  <>1  a slati.sticallv  iiomoni  lu'nus  material  mav  ii"i 
lioKI  li.r  larlain  luhamnr.  I lew  e\ar,  when  it  (lot's  apply  and  when  the  piineiph 
ot  local  action  apiilies,  we  can  use  eorrelat  i\ e models  supported  l)\  adecpiale 
testinu  till'  stiaietural  meehanies.  ( iperatioiially,  this  inav  I'e'juire  one  or  nioi  e 
pi(.'  processors  to  obtain  a eorrelatite  model  and  one  or  more  post  processors 
to  obtain  results  at  the  constituent  le\el.  (liven  this  premis(.  , the  iirisent  etfort 
loeus(.'d  on  tk'\ elopin;;  a eompidalioiial  system  lor  all  those  correlative  motlels 
that  can  b(.'  i'haraet(.'ri .'cd  by  state  variables.  Ihis  is  a very  lar^c  class  and  will 
allow  the  numerical  analysis  ol  many  composite  structures  whose  inelastic 
behavior  differs  from  nu'tal  plasticity.  I'lie  approach  also  allows  the  model  to 
function  as  its  own  pre-processor  and  synthesi/e  eorrelati\i'  models  Irom 
I'onstiluent  inela.stie  propertii,'S  and  the  bi'liax  ior  at  interlaces. 

The  magnitude  ol  the  computational  pniblem  associated  with  inelastic 

Ihri'c-dimensional  analyses  a))proaehes  that  ot  (.■haraet<.‘ri/ inp;  tln'  materials 

Ithvsies.  Several  (.xei'llent  suiar'V  paiiers  are  available  as  well  as  ease  histories 

12 

of  recent  applications.  The  computational  ajiproaeh  taken  in  thi'  |)resent  |)aper 
is  based  on  the  luemise  that  better  computers  are  more  likely  than  bi'tter  numerical 
methods.  C'onsr'ciuently , a computational  meth(Hl  was  adofJted  that  takes  advantage 
of  'Ahat  most  scientific.'  computers  rio  Ix'St,  vector  iiroeessinj;.  It  is  an  adaption 
of  thi'  eoniunate  gradient  algorithm  that  avoids  eoimeetivity  optimi/ation,  that 
uses  relatively  little  eori'  and  w!v.;se  cycle  time  grows  at  most  linearly.  One  of 
the  method's  disadvantages  in  linear  problems,  it  is  inefficient  for  multiple  load 
conditions,  can  become  an  advantage  in  nonlinear  problems  when  the  iteratif'n 
for  the  next  load  increment  is  started  from  the  solution  for  the  previous  increment. 
A second  as))ect  of  the  computational  problem  is  that  of  data  generation  for  the 
discrete  mathematical  model.  This  problem  ree(.'ived  considerable  attention  in 
the  original  development  of  PA  I'C' I IKS- III  ami  resulted  in  a parametric  cubic 
tiKHk'ling  system  based  on  construction-in-eontext . The  pi  esent  effort  addeil 
new  geonu'trv  constnu'tion  operations  to  the  system.  de\c'b'ped  the  paramelri- 
/ation  necessarv  for  modeling  strain  singularities,  di'vi'loped  intra-element 


k 


V)roiKTty  muiic'lin^r  aiul  developed  a eommon  ini)ut  format  for  all  nonlinear  material 
))ropertic‘.s.  Applications  of  these  features  are  made  to  a carhon-carb(jn  unit  cell, 

to  an  interior  crack  problem  and  to  the  triaxial  response  of  a graphite  bar  tested 

. 2 
by  Jortni’r. 
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MA  ri;i<iAi,  Mt  »i)i'  1.S 


. 1 .SlaU'  \'ai  i:il)U'  MihIcIs 

rhc  i\ (■  is  a ^i  iu  i'al  i-ompulatioiial  mndi  l lliaf  accnunts  lor  the 

anomalous  incl.istie  hi  l>a\ioi-  ol  eomposite  malci’ials,  'I'alilc  1.  The  Ini.sic-  assump- 
lioiis  arc  that  tlie  material  is  slatistieally  homo>>eneini.s  and  that  at  the  mac-roseopie 
level  inelastii-  hehavior  is  stric  tly  l<»eal.  It  is  also  tacitly  assumed  that  this 
behavior  is  dederminist ie  which  may  be  the  weakest  link  in  the  chain  ot  assumptions. 
Physically  based  statistical  models*  * have  beam  very  sueeessl'ul  in  explaining  the 
extensional  beha\  ior  ol  graphites  associated  w ith  tJie  nueleation  and  gi-owth  ol 
mieroeraeks.  However,  additional  work  is  reiiuired  to  inedude  a shear  mechanism 
and  to  rc'duee  their  computational  expense  before  tliey  can  serve  as  the  basis  for 
a general  eom|)utat ional  model.  Mathematically  basi-d  statistical  tiu'orii's  have 
not  been  nearlv  as  sueeessful. 


Table  1 .Anomalous  inelastic  behvior  obseixed  in  eom|5osites 


Behavior 

.Metals 

Composite'S 

Yielding 

Independi'iit  of  hydro- 
static st  ress 

Can  bi'  di'pendi'iit  on  hydro- 
static' .strc'ss 

Inelastic 

Mechanism 

C rystallographic  slip 
Twinning 

.M  ic  roe  rack  ing , in te  r st  itia  1 
slip,  mierobuckling 

Inelastic 

.Strain 

Constant  volume 
plastic  strains 

Inelastic  strains  c;m  cause 
volume'  change's 

Unloading 

( K’curs  at  elastic 
modulus 

Can  occur  at  secant  modulus 

Failure 

Convex  failure  sur- 
faces 

Concave  failure  surfaces 
can  occur 

The  approach  taken  a.ssunu>s  simply  that  the  stress-strain  law  is  some- 
function  of  one  or  more  stress  or  strain  state  \-ariabk‘s.  This  assumption  can  be 
list'd  for  both  incremental  and  total  stress-strain  laws  but  tlie  pri'sent  effort  is 
limited  to  the  latter  and  to  three  state  variables 


- 4 - 


V,  I- . , , (\', , V ) a cr,  , (.'{| 

2 1 ikl  \ I 2'  II  kl 

'I'lii'iv  is  :i  simil:iril\  li>  tiu’  ni>nliiU':ir  stri’ss-strnin  mudcl  ol  N’ucbiM'  , l)Lit 

liis  mi  Ilk  ■ I |•l•';|  i-ifts  llu-  Inrm  n|  I-  i|u;iti"n  i 1 1 :iml  ik  n's  ;i  Ik  i\\  tlir  ( • n'M  iciints  in 
K(|ualii)iis  i2i  and  i:!i  in  he  indrpi'ndi'iil  l"iK'linns  n|  llu-  slati'  ol  tlu-  mati  I'lak.  I'his 
is  lU'i'i  ssa rv  to  ai-ooun(  lor  ohaiUfii  s in  tin  do^riu'  ol  anisotropy  during  iiu-lastic 
straining.  * ^ \ great  manv  i-orn-laliM'  inodi-ls  including  Hatdort,  loncs-Nclson, 

Weilcr,  llatm-Tsai  ;uul  the  classical  deformation  models  can  be  represented  using 
K(|uation  ill.  It  allowed  the  design  ol  a nonlinear  mati-rial  module  eiu-ompassing  a 


\arietv  of  lormulations  all  using  one  input  tormat.  All  moduli  and  stati.-  tunctions 
are  inleriiolated  using  piecewise  cubies  which  are  callable  ol  approximating  any 
continuous  funclion  as  closely  as  required.  This  a\oids  the  interpolation  probk-ms 
associateil  with  exponential  approximations  in  the  lones-.\elson  model,  rlu-  stati' 


\ariable  used  in  this  model  is  the  strain  energy  densdv,  W,  so  that  the  F..,  , 

2 i|kl 

are  the  ('  , , and  llu-  nonlinear  eiiuation 
1 jUl 

W ('  (VVif  * f | (loncs-Nclson)  ill 

lllsl  II  kl 


must  be  soh’ed  for  each  strain  stale.  I'hc  Hatdorl  model  is  a strain-stress  lormula'ion 
that  is  \ery  simple  and  accurate  for  transi I'rselv  isotropic  graphites. 


It  usi  an 


r ;i  state.'  variahli'  llial  e aii  lie  com  iTti'd  to  a sli'aiii 


I 

I . I 

1 


I'fli'i'liNt'  stri-'ss  paraini'liT  loi 
tiH'inulatii'ii  l>y  solving  a nniiliiu'ar  I'emation. 

__  *) 

a “ 1 , , a.  <J,  , ( Maldnri  > 

iikl  i|  Isl 


t'.  '<^•^11  , 

ilinn  klsl 


M 

inn 


,\1 


(.11 


riu'  moduli  ai'o  linearly  ilc  pi  ndcnl  on  o in  this  modi'l.  in  order  to  aeeount  toi' 
tensioii-eompressioii  heha^  lor  tliilerenees,  llaldoi't  uses  an  ad  hoe  proeedure  based 
on  llie  sit’.n  <’t  individual  stress  eoniponenls.  t'niortunately,  this  pi’oeedure, 
when  used  with  I'.'qaation  olteii  will  k'ad  to  oseillafiiif;  die  i'rf4i'ni-i'  between  two 
stress  slates  havini;  tiu'  sanu' tr  but  dilt'erent  si^ns.  To  aeoid  this,  the  si^>ns  are 
held  eonstant  elui  int;  eaeh  inaior  cyc  le  in  PAT('1II!S-1I1.  These  ditlieulties  are 
rewardeil  bv  surprisingly  aeeurate  a^;reement  with  the  extensice  biaxial  and  Iriaxial 

9 

data  obtained  by  Joriner.”  Numerical  results  obtained  usinn  i jjualion  (.ii  are  com- 
pared with  his  data  later  in  the  paper  for  tria.xial  loading  at  constant  strain  ratios. 


The  Weiler  model  also  uses  an  effeetice  stress  i)arameter  for  a state 

varitible  but  allows  the  T , , to  be  dependent  on  the  efteetive  stress  to  account  for 

i |kl 


changes 


-2 


cr 


IcJ  ' C7  a,  1 I VVeileP 
n k 1 


(b) 


in  anisotropy  with  strain.  The  in-ocedure  used  by  Weiler  and  others  is  based  on 

constraining  the  effeelive  plastic  work  to  be  epual  to  the  uniaxial  plastic  work  in 

each  component.  Aj)pl ieations  of  these  mot'els  ha\e  used  functions  that  make 

the  effective  sti'ess  insensitive  to  hydrostatic  t)ressure,  but  this  is  not  a reciuiremenl. 

Hvbicki''^  included  plastic  volume  change  in  an  V.  model  for  a carbon  material. 

1 ) 1 

•JTA,  several  years  ago  and  obtained  good  correlation  with  data  for  a ))ressuri/.ed 
test  cylinder  under  axial  load,  lie,  too,  used  an  ad  hoe  proeedure  to  account  foi- 
tension  compression  differences. 


The  correlative  models  described  have  been  used  principally  lor  partic- 
ulate eom()osiles.  It  is  also  i)ossible  to  represent  most  eori'i'Iati\ e models  used 


J 


for  fibor  roinfori’od  oompositos  with  Kquation  (1).  Consider  for  oxami)lc  llio 
lialm-'I'sai  modcd  in  whioh  tlu>  nonlinoai-ity  lor  a laminar  is  shcxu'  dominated. 
In  tliis  i-ase 


a , 

1212  12 


(Hahn- Tsai) 


and  the  cxtonsional  f,  , , an-  indc'poiuk'nt  of  V., . riiero  are  essential  clilTerenees 
ijkl 

between  this  modi'l  and  those  usin';  strain  energy  for  although  at  small  strain 

“ 21 

levels  both  tit  the  data  lor  graphite-epoxy  laminates  reasonaljly  well. 

2 . 2 Spatial  X’arialile  Models 

The  objei-ti\e  is  to  model  property  distributions  as  aecurately  as  defor- 
mations so  tJiat  modi'ling  deUii!  is  eontrolled  l)v  representatitm  of  the  output  rather 
tlian  the  input.  This  reiiuinmient  is  another  instaiu'e  in  whieli  eomposite  material 
modeling  differs  from  metals  both  for  mierostructures  and  struetures.  carbon- 
earbon  unit  eell,  for  example,  ean  l)e  modeled  with  one  variable  property  element 
or  eight  eonstant  property  elements.  A helically  wounti  structure  with  changing 
radii  has  properties  that  are  continuously  changing  and  is  also  most  efficiently 
modeled  with  variable  property  elements.  Ttie  approach  tal^cn  in  PATCIIES-III 
is  to  use  the  same  parametric  cubic  functions  for  l)oth  geometry  and  physical  data. 

It  is  a generalization  of  the  isoparametric  approach  ami  requires  the  development 
of  new  modeling  tcchniciues  for  properties. 
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Hecently  Henshell  demonstrated  how  variable  mesh  point  spacing 
(parametri/ation)  can  be  used  to  greatly  improve  isoparametric  models  near  a 
strain  singularity.  The  same  technique,  changing  the  parametrization  of  the  geometry 
moilel,  can  be  used  to  model  the  step  in  material  properties  at  an  interior  bimaterial 
interface.  In  this  instiuice,  TTgure  1,  we  induce  an  inflection  point  in  the  property 
model  at  the  interface  by  eonstructing  a geometry  model  with  :ui  inflection  point  at 
the  interface.  Consider  a coordinate  function  '/(^)  in  algebraic  format 


;uul  Uk‘  ('‘Hulitions  uivi  iullfclion  t onstruint 


/((») 

/.l,  s,  s, 

/([,  0 :i  Sj  ^ • 2 sj  S,J 

/(i)  u = <;  • 2 s.,  f'Ji 

wluTc  ttu'  step  is  ;it  V.  /(i).  Next  solve  l•:()u;lti^>l)s  ('.M  in  terms  ot  the  as  yet  unJ-;tU)wn 
5 where  A'/  ' '/(M  - /.(<')• 

s.j  A/  (:ir  - i) 
s.,  - -:!A/(  far  - af  ■ i) 
s.j  - aA/r  ' (a?"  - af  • i) 

S = /.(O)  (!••) 

I 

riu-  location  oi  llu-  step  or  interlace  in  parametric  space,  r can  now  be  found  by 
sohing'  thi'  cul)ic  equation 

- a (A^^AX)  a (Ax/Ax;  ? - (Ax/Axi  = o /iD 

for  the  I’f'ot  in  the  interval  a ^ ^ < I wliere  AX  = X - X(0)  and  Descartes’  rule  ot  signs 
guaranti'cs  theme  will  be  a root  in  the  inters  al.  Finally . a parametric  cubic  for  a 
property  component  I’  (^)  is  oljtaini'd  Ijy  imposing  the  boundary  eonditions  P ((')  P^^, 

P (1)  Pj  and  P di)  P (1).  Other  choices  are  possible,  for  example,  the  value  of 
P (^)  could  1k'  prescribed  at  four  points,  Init  this  will  cause'  P t^i  to  fall  outside  the 
inte'rval  [l’  . P over  portions  of  the  interval  0<  1.  This  may  be  necessary  if 

i '>  . . 

the  area  under  the  curve  is  a key  parameter.  One-  might  keep  the  P (01  = P (It  = 0 
e'eeiistraint  and  choeise  the  value  of  1^  (0)  such  that  the'  intt'gral  of  P (^)  on  the  inte'rval 
()<f<f  rives  tlie  correct  area,  P AX,  with  the  cen-responding  choice  feu'  P (1).  In 
ge'iieral,  the  eU'term ination  of  the  parametri/.ation  and  the'  intra-eh'ment  prope'rty 
moelel  form  a nonlinear  pre)gramming  prol)lem  otu'c  an  accuracy  criteria  such  as 
le'ast  seiuares  is  e'stablislu'd.  Ilovve'ver,  a great  de'al  can  be-  accevmplislu'd  with 


simpli.'  aiKilysc's  ol  tlu‘  tyiu-  prosi-nk'd.  Paraimdrizations  lor  a strain  siunularity 
at  Z (0),  tor  i-xamplo,  can  1k‘  obtained  without  solving;  a cubic  and  sc\’cral  are- 
listed  in  Table  '1. 

Tal)le  - Strain  singularity  para  met  ri/.ations 


Z(^) 

ZiO) 

Z(1  ;j) 

Z(2/3| 

Z(l) 

• ) 

0 

L/‘J 

iiyit 

L 

k:i 

s 

0 

L/27 

- 1.-27 

L 

s 

0 

■IL/27 

20  L/27 

L 

C'onsicU-r  ni-\t  the  modelinf-  of  a fiber  reinforced  structure  in  which  the 

orientation  of  the  mati  rial  axes  chan)4i-s  continuously  \\ith  respect  to  a reference 

frame.  An  intc-restiiif-  i-xanipU-  is  rosette  construction  used  in  rocket  nozzle 

structures  in  which  the  lilu-rs  spiral  about  an  axis  with  changing  radii,  Figure  2. 

2;t 

I’agano  provides  a ih-tailed  account  ol  the  relations  between  coordinate  tvames 
and  shows  that  while  axisymnu-tric,  the  jiropc-rties  \ary  in  the  radial  dirt-ction. 
Also,  all  21  elastic  constants  are  non-zero  for  this  mati-rial  in  cylindrical  or 
rectangular  coordinates  making  t-fficient  property  modeling  particularly  important. 
The  transformation  from  mate-rial  coordinates  to  cylindrical  coordinates  in  this 
case  is  ck-termined  by  a rotation  constant  followed  by  a rotation  a a (r)  where 
r sin  rt  constant.  The  transformation  to  re-ctangular  coordinates  simply  rec|uires 
ailding  a rotation  0 to  a.  Intra-clement  prope-rty  modeling  in  this  case-  requires 
only  the  ply  prope-rties  anil  the  spatial  distribution  of  the  Kule-r  angles.  The  ori- 
ginal I’.ATC’ll  FS-lII  syste-m  was  eie-signed  for  such  input  but  did  not  anticip;ite  the 
lu-e-fl  for  vai-iable-  Fuler  angle  elata.  As  a i-esult,  the-  spatial  \ariation  ot  all  21 
elastic  const-ants  had  to  be-  nH'deleel  individually  using  a pre-proce'-ssor  for  a one- 
ele-me-nt  modi-1  of  a rosette-  cylinde-r.  The-  same-  moeiel  coulel  ha\e  been  createil 
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Irom  ttu'  sp;iti;il  variation  oi  threo  Kulor  anulos,  two  ol  whii-li  aiv  constant.  Com 
parisons  Ik'Iwooh  tlu‘  paramotvic  (.ailjii’  propfi'tios  ami  tlio  analytic  sliowcd  dillcrcnccs 
of  K'ss  than  one  i)crccnt.  Coini)lctc  results  from  this  study  will  lie  |)resenti-d  in 
a later  paper.  Finally,  it  shouhl  he  reinemhered  that  variable  propi-rtv  modeling 
adds  no  new  ik‘^;rees-of- 1 reeilom  to  tiu'  anal\  sis  model  and  can  reduce  flu'  number 
of  idements  ri-quired  for  some  composite  imderials. 


1 


0 


C'(  )M  IM  rA  TK  )XA  I , M(  )1)K  LS 


1 (ieoincUy 

I'hf  basic  constructions  and  [iroiK-rties  foi-  ixiramctric  cubic  line,  surface 

21  , ].■{ 

and  volume  models  may  be  tound  in  the  work  ot  Coons  and  others  . 1 hey  pro- 

vide a data  base  foi  compuli'i  aided  p;eometric  ilesiy,!!  comprehensive  enoujth  to  be 

2f> 

used  for  pai'ts  definition  in  manufacturing;  and  they  now  are  used  extensively  in 
industry  here  and  abroad  to  <lefine  external  suidaces.  Their  reliability  tierives 
from  modeliii};  in  parametric  space  where  there  are  no  proljlcms  \Mth  asymptotic- 
slopes  and  their  accuracy  derives  from  using  llermite  polvnoinitils  which  interpolate 
a function  and  its  derivatives  with  the  smallest  possible  error.  The  price  for  these 
qualities  is  twelve  coefficients  for  a line,  fort\ -eight  lor  a sui’face  patch  and  one 
hundred  ninety-tw'o  for  n colume  hyperpatch.  The  dalti  generation  pr"b!em  for  these 
coefficients  can  be  solved  using  the*  construction- in-context  a iproach  developed  ori- 

gimdly  tor  I'.ATCllES-IIl.  nils  system  uses  a variety  ot  TI.M-. . l’.-\TCli and  HP 

direc-tives  with  cross  referencing  to  c-onstruct  the  geoiiu-tr\  model.  Two  con- 
struction opertilions,  ruled  volume  and  outline  surlace,  are  illustrated  in  Figure  d 
and  a listing  of  the  currenllv  available  options  is  provided  in  Table  d . The  outline 
surtace  construction  opertUion  shown  in  Figure  d consists  in  having  the  computer 
move  an  outline  curve  along  a base  cune  with  a fixed  orientation  relative  to  either 
a global  Cartesian  frame  or  the  local  Frenet  frame  of  the  base  curve.  In  developing 
this  option,  it  iH-came  ob\  ions  that  in  many  instances  it  would  be  desirable  to  changf' 
the  initial  orientation  of  the  outline  curve.  This  feature  was  provided  by  allov  ing 
an  inilitil  tran.'4'ormaiion  that  leaves  the  file  copy  of  the  outline  curve  unchanged. 

In  retro.stiect,  manv  o|  the  original  PATCHFS-III  directives  would  have  benefited 
from  such  a feature  which  functions  as  a modifier  in  the  language  implicit  to  the 
present  approach.  The  benefits  of  having  even  a primitive  language  have  been 
substantial  for  data  generation.  The  system  functions  as  its  own  pre-processor 
and  allows  shapes,  like  those  shown  in  Figure  4,  to  be  created  using  roughly  ten 
input  directives  (cards)  per  model.  However,  the  models  created  contain  sevi’iiil 
hundred  coefficients  that  completely  describe  the  geometry  of  the  figure.  The 
Frenet  frame  fora  line  /(^t,  for  example,  can  be  computed  directly  from  the 
parametric  cubic  mixlel 


- 11 


(Ili  I 
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w !km\'  tlie  (lot  iiiilii-;iu-s  (liHi'icntialinii  with  resport  to,  the  paraim-lr'U'  (•(.oi'dinaic. 
It  is  also  a |■olltinl'  matti-r  to  eompiiti'  surface  normals,  areas,  volumes,  eur\atuii-s 
anil  in  sliort  anv  neometrie  iimpertv. 


I'ahle  ;i  (ieomelrv  eonst  ruetion  lii  rei  tn  es 


1, ARC  PC 
UNFP 
1.  IN  lies 
1,1  Nil  PC 

PA  rciiR 

P.ATCIb  di 

PA'l'Clll. 

PATCIIO 

PATCIKQ 

I'AI'CIIR 

111  'H 

llPl, 

'IPllllX 

IIPN 

IIPP 

IIPR 

11P2PAT 
lIpc.PA  r 

SCALP 

SCALPll 

lMn\  11 

‘ Mnemonic  Suffixes: 

A .Algebraic  , 

R Ceometric 

CS 

Cubic  Spline,  1,  Line, 

N ■ Normal 

P Point 

PC 

JMrametric  Cubic  , 

(Quadrilateral  , II  RoLation  , 

: 1 . 2 Fi n ite  F lenient 

The  decision  to  use  tlu'  sixty-four  point  isoparametric  finite  element  for 
the  present  study  rests  on  its  modeling  efficiency  for  anisotropic  materials,  for 
palhologieal  shapes  and  for  general  boundary  conditions.  There  are  many  cases  in 
which  these  (lualities  are  not  required;  however,  the  focus  here  is  on  the  anomalous 
behavior  of  composites  particularly  in  highly  stressed  states.  In  this  rcgtird  it  is 
important  to  reali/e  that  material  anisotropy  affects  matrix  conilitioning  as  strongly 
as  geometric  asiject  ratio.  Consider,  for  example,  a heated  disk  with  isotropic 
properties  (11/ C 2.<il  and  with  highly  anisotropic  proi^erties  (H/G  17.2a)  in  which 
the  temperature  varies  (luadralically  with  I'adius.  This  change  in  material  anisotropv 
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i-aiisoil  tlie  iuim))er  ol  clometUs  required  lor  ^;u(k1  stresses  to  double  as  I able  -1 
demonstrates.  The  lar^e  hoop  stress  ratio  at  r 2H/3  is  near  a 0^  0 point 

Table  4 Anistotropv  elTeets  on  stress  aeeuraiy 


Isotropic  .Anisotropic  Anisotropic 

Two  Ideinents  Two  Kle meats  Four  Klements 
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1 . 

M2  1. 
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•) 
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:dJ7 
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• 
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• < 

1 . 0D‘J 

. 

Hatii 

1 of  Cl 

imputed 

1 Stress 

to 

exact  stress 

• Th.' 

c.xaci 

St  I’CSS 

is  zero 

at 

this 

point 

and  IS  simpl\  an  mdicatioi  of  larne  differences  in  small  numbers.  Similar 
studies  on  the  effect  ol  hiuh  a.spect  ratios  and  degenerate  element  shapes  indicate 
the  sixtv-four  point  element  has  a broad  band  over  which  it  can  maintain  good 
stress  accnracv.  In  order  to  realize  these  benefits  many  practical  problems 
associatc'd  with  large  element  matrices  had  to  be  solved  in  the  original  PATCHES-ni 
system  and  the  added  Inirden  of  material  nonlinearity  required  similar  efforts. 


An  t‘\:imiii:ili"ii  <>r  tlu‘  |•(■la(iv^■  costs  of  ^cncratin^  an  i-lcnuait  stilliU'ss 
matrix  \s.  ^ciu'ratin;.;  an  rlt'ini'nt  Ihcnnal  lo;ii!  vector  lui-  ;i  carljon-carhon  material 
raiiKc  Irom  10:1  to  over  ‘JO:  I ile))<m(lim4  on  the  "ciimetry.  'Hiese  I'atios  ari'  \'ery 
liis^ii  in  comparison  to  lowao’  or»ler  elements  anti  stron^lv  su^jttxst  the  use  <j1  a f>st*utlt)- 
loree  I'orimilation  for  mali-rial  nonlinearity.  When  ilireel  matrix  solution  mi'thod.s 
are  usetl,  another  factor  in  favor  of  this  a!>nroaeh  is  the  ratio  'if  matrix  assemMy 
and  deei)in|)osition  (ime  to  tlu‘  forwaril  liael.ward  solution  time  which  is  t\piea!ly 
(■>;!.  The  eompul at iona  1 factors  in  lavor  <'1  :m  initial  stress  formulation  are  laster 
I'oin  erv't'iiee  and  vrealer  stability  at  hi;;!)  sti'ain  It  vels.  Stutlies  of  the  relati\e 
eonvervimee  rates  by  ilavner’’  int!ii,'ate  |•||uvhiv  If  • 1 ratios  in  favor  tif  (he  initial 
stress  me-thotl  \\hieh  is  nt>t  nearl>  t laiimh  to  mat-f  the  methoi!  ettmpetitivc'.  Staijilitx’ 
\^il]  not  be  a faettir  for  snnill  stiains  anti  a peritulie  initial  stress  cycle  can  always 
bi  tahen  if  neeess.i rv. 

Tht're  are  many  tlerivations  of  the  matrix  ctjuatitins  assticititetl  with 

the  pseutio-foree  nudhotl  ami  they  nceil  not  be  i-epeateil  here.  The  matrix  e(|uatitin.s 

2(1 

usetl  in  the  present  sfiuly  are  very  similar  to  tlmsc  uscal  by  Ilavner. 

|K(V,  T)|U  F 

IKd/niU  r 1K((  b T)  - Ki\', 'I’ll  r dSi 

which  leatls  tfi  the  reemrsion  relatitnis 

|K(i),Til  r , F IKtO.Tl  - K l\  ,Tlj  F 

' ’ ~n  I ~ n ~n 

|Kdi,T,|F  F ‘i  (111 

~n  1 ~ ~n 

The  pseudo-force  term  is  tibttiined  by  inte;;r-atiii{i  the  strains  frtun  cycle  n aj>,ainst 
the  difft'retu f in  propertii's  relerrini;  to  Fijiure  i 

b'  f ^Ht  ^A('  1 ( tl\'  il  'd 

~n  /•  n ~n 

\ 

where  (H)  transftirms  mesh  ptunt  tlis]ilaeemcnts  F into  strains. 


M - 


A SLiltii'k’at  L-omlition  loi-  (.•<)n\frp;c-iR:c  oi  tiu'  mcthoii  interred  from  IUiviut's  aiudysis 
is  tliat  l')|  must  hi'  monotoiiii'  in  \’.  It  should  In-  nuti'd  that  in  material  modids 

V.  ith  iiu'laslii-  \ oIumi'  i'liaai;i‘s,  tlif  hody  torccs  tor  tluTinal  stri'ss  prohloms  will 
also  I'haiigi'.  I'hi'si'  chairai'  oan  he  accounti'd  tor  in  by  using  tlu*  mechanical 
strains  in  taiuation  (lai.  I'hc  limited  cxpi'ricnci'  to  date  with  the  method  has  enjoyed 
rapid  eomergimee.  As  more  complex  prolilems  are  attempted,  it  may  be  ni'cessary 
to  utilize  more  etticiimt  ri'eursion  relations.  Another  method  that  was  considered 
in  the  study  easts  the  prol)lem  in  the  term  ot  a tirst  oi'der  ditterential  ecgiation 
using  the  residual  lector  tormed  trom  tiiiuation  (lb. 

K -|K|h  ■ K g (Ki) 

u . H - g a 7 ) 

Tlu'  dot  indicates  a time  deri'.  ative  w here  time  increments  are  synonymous  with 

load  Increments  and  the  sealai  X is  equivalent  to  an  over- relaxation  laetor. 

27 

Stricklin,  et  al.  use  this  approach  to  derive  a sell-correcting  procedure  that 
worked  very  well  in  tlieir  nonlinear  material  api)lications. 

3 . 3 .Mat rLx  Solution 

Thi'  third  major  eomimtational  problem  after  data  generation  and  matrix 
generation  is  matrix  solution.  The  use  of  a (i-l-point  finite  element  leads  to  matrix 
equations  that  are  relatively  dense.  It  is  not  unusual  for  densities  of  33  percent  to 
occur  in  matrices  well  over  dimension  looo.  The  original  linear  code  development 
concentrated  on  the  first  two  computational  problems  and  used  existing  (,\.ASTKAN) 
matrix  routines  to  solve  tlie  matrix  eriuations.  As  a result,  tin-  system  is  efficient 
for  probl(>ms  up  to  alji'ut  rlimension  IdOO  and  looses  efficiency  for  larger  matrices 
until  eventually  saturating  the  computer.  The  need  to  solve  similar  dimension 
nonlinear  matrix  eijuations  led  to  a change  to  an  iterative  solution  method.  The 
advantages  of  this  change  are  1)  elimination  of  large  matrix  files,  2)  freedom  from 
connectivity  optimization  problems,  3)  relatively  small  core  requirements  for  large 
problems  and  t)  heavy  vector  processing.  The  disadvantages  are  possible  shnv 
convergence  and, in  lineai'  problems, N multiple-load  conditions  are  N times  as 

•> 

expensive  as  one.  The  iterative  method  selected  is  the  conjugate  gradient  algorithm" 
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;iii(l  ill  |irnpci-ly  scaU'il  cooi-ilinaLc.s , it  (.•nnviTpcs  in  iar  IcwtT  > yclcs  than  thr  ilimcnsion 
oi  tlu‘  inalftx.  Uc  l(MTin:4  to  i:(|Uatii>n  (1*m  ti>r  the  ri'sidual  It 


I'.,  i!'~ 

I r I I> 
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I'  Iter  ,1  • i'") 
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wlu'ia'  i1k‘  sfah-rs  t,  anil  ft.  arc 
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t I U.l  “ / l‘.  |K1  l>. 

1 ~i 

il!i  r 'li,  ' 


dill 


lliiH  is  liu’  linear  iorin  ol  the  al^erithin  aiipri'iiriatc’  to  the  si  hit  ion  o!  I'oualion  (1  1). 


It  is  also  possible  to  usi’  thi‘  nonliia-ar  \a  rsion  ol  tin  al^toritlim  to  sol\e  laiuation 

T 

(K!)  liireetly.  In  tliis  I'ase,  t.  is  the  smallest  positive  root  ol  < d J: 

however,  past  exiierieiiee  lavors  sueeessivt'  elastie  solutions  and  this  was  done. 
Iteecut  applieations  of  the  lim  ar  alt;orithin  to  HI)  composite  material  proldems  by 
Dana”'*  always  eomerai’d  to  four  places  in  less  lhan  N'/o  cycles  lor  problems  ot 
dimension  N obO  to  over  2non.  Applications  in  the  present  study  hive  b-'en  to 
small  onc'  and  two  elc-menl  models  recpiirine  X/2  cycles  tor  similar  aceiiracv. 


Idficieiiey  is  also  a function  of  the  c’osi  jn'r  cycle  aiul  this  is  whc-rc  the 
vector  procersinp:  efficiency  ot  scientific  computers  1rI])s  iterative  alttoritlims. 
l-lven  till'  CDC  (id'M)  can  lie  made  to  compule  dot  products  very  elficiently  t)y  tahina 
advantastc  of  multi))le  arithmetic  units,  rnfortunately , older  computers  like  the 
U.MVAr  110s  do  not  have  tins  feature  and  Ijecause  of  tlieir  short  :i2  tiit  word,  all 
arithmetic  must  lie  vionc  in  double  precision.  Several  attempts  to  use  mixevl  mode 
arithmetic  on  a L'NJ\  A(’  I IDs  were  unsuccessful.  T!ic  operation  most  sensitive 
to  round-off  error  is  the  transtormntion  of  an  element  stiffness  matrix  from  neo- 
metric  lormat  to  point  lor  mat.  The  reason  is  the  lai'ne  diflereiice  in  (he  mannitude 
of  disidaeemeiits  and  displacement  derivatives.  !•  Ortunatidy  this  can  lie  avoided 
by  transforminn  directly  Irom  aln“braie  format  t"  point  format,  Anotlver 
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faclor  atTi'ctiiifi  ryi'K'  citicu-ncy  is  the  amount  of  data  transfer.  In  the  pri'sent 
apin-oaeh,  only  the  redueed  element  nuitriees  are  transferred  into  eore  eaeh  eycle. 
older  vi-rsions  of  tlie  pi-oj;ram  transferred  the  redueed  struetural  matrLx  into  core 
eaeh  eyele  and  this  was  more  expensive  primarily  because  tlie  matrix  had  to  be 
unpaeked  a few  eolumns  at  a time.  As  a fringe  benefit  <jf  using  element  matrices, 
tlie  data  transfer  per  cycle  grows  linearly  with  the  number  of  elements  in  the  model. 
The  overall  eyele  eosl  then  grows  at  most  linearly  sinee  vector  processing  also 
increases  at  most  linearly  witli  the  number  of  elements. 

A constant  issue  in  itc’rative  metliods  is  the  convergence  criteria.  Tlie 
first  parameter  to  eonverge  is  the  energy  with  the  maximum  modulus  displacement 
component  usually  a close  seeoiul.  When  these  two  parameters  have  converged  to 
seven  places,  the  stresses  and  strains  have  about  three-plaee  accuracy.  Several 
adtlitional  cycles  are  reciuired  to  produce  stresses  luid  strains  that  agrc-e  with  the 
direct  solution  to  six  places  (see  figure  <b  :uid  this  raises  the  issue.  Should  these 
additional  cycles  lie  executed  to  increase  the  accuracy  of  the  least  significant  stresses 
and  strains  and  provide  eonsistcncy  with  direct  solutions?  In  linear  problems  this 
is  done  because  the  increase  in  cost  is  modest.  However,  in  nonlinear  problems, 
this  expense  is  harder  to  justify  since  the  internu-diate  solutions  have  no  function 
other  than  t<>  provide  accurate  state  variables,  like  strain  energy,  so  that  accurate 
material  properties  can  be  computi-d. 
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1.  I) 

1,1  ( '.1  rhi)i)-t’:i  j 1)011  l iiit  Cell 

riio  prooi'i-siii'r;  ot  'il)  c-aibon-c-:irl)on  malc-nal^  idoullv  results;  in  an  orllio- 
>;onal  ai  i-a\  oT  niun's  in  whii’ti  the  opL-n  r(.‘nionh  art-  lillod  v.iUi  a I’arbon  inatnx  inalu- 
I'ial  (lui’in^  diMisilicatinn.  fno  i (,‘pt‘li tivo  voiunu*  eltMiienl  I'of  die  nialerial  lias  lliree 
planes  o(  svnunelnv  and  I'’ti4tire  7 shows  selieniaticallv  one  oetanl  ol  the  so-ealled 
iinil  ('i‘ll.  .\hK’ fo.'ieopi [)roperlies  lor  this  material  arc*  usually  liasc^d  on  an  aiuilysis 

ot  the  unit  eell  with  uniform  traction  or  dis|jla<  cmeiit  hoimdarv  coniiitions  as  des- 
cribed by  boss'*'*.  A constant  property  finile  element  model  of  this  coinpusiU*  would 
I'eciuire  a minimum  of  ei};ht  elements  and  the  rc*fercnccd  analysis  used  twenty-seven 
elements.  This  same  unit  c*el!  wa.s  modc'led  with  only  one  cairiable  properly  element 
usinp,'the  techniciues  de.scribcd  I'arlier.  I’he  consUluent  properties  lor  the  libei 
buiklles  and  matrix  material  arc  p,iyen  in  Table  a and  a more  coiu|)Iete  descriptiijii 
of  the  composite  is  aivcii  by  Koss‘*^\  The  results  from  the  one  element  model  are 
in  remarkably  u;ood  a>;reement  with  the  results  Irom  other  tinalyses  presented  in 
Table  I).  There  are  several  commc*nts  that  need  to  be  made  about  these  results  to 
maintain  perspective;  ill  the  orthogonal  nature  of  the  weave  allows  the  macroscopic 
(jroperties  to  be  computed  from  coarse  models  (2i  the  comjHiter  plot,  I’igure  s,  ol 
the  free  thermal  expansion  of  the  one  element  l’.\'TCIlKS-in  model  shows  mesh  lines 
at  f L/'ii  to  enhance  viewiiit;  and  (d)  i mperl'ections  in  the  unit  cell  can  result  in 
large  differences  between  ideal  and  real  malerial  properties. 

(liven  this  perspective,  one  conclusion  to  be  dr.awn  is  that  without  varial)le 
property  minleling  an  analysis  using  the  pre.sent  element  would  cost  over  an  order  of 
magnitude*  more  ;md  not  substantially  change  the  mechanical  pmperties.  A second 
interesting  conclusion  is  related  to  the  poor  estimate  of  a 1 1 and  the  sensiti\ity  ol 
the  macroscopic  properties  to  constitueni  properties.  The  parametric  cubic  property 
modeling  used  l'(ol  1'^,.  P(l)  P,  and  i^(ol  P(l)  - 0 boundary  conditions  which 
tend  to  lower  the  high  modulus  constituent  data  and  raise  the  low  modulus  data. 

A sensitivity  study  revealeil  that  all  was  primarily  a function  of  the  longitudintil 
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Table  f)  Tnil  cell  dimensions  and  constituent  properties 
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reference  30. 

d'ahle  I)  Unit  cell  macroscopic  property  comparisons 
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•SAP  results  from  reference  30. 


a uT  tlic  '/^  iiht*r  wliich  is  viTV  small.  In  this  ctise  the  p;tfametrir  culjif  modol 
cltortivolv  liouhlml  the  a ^ of  tho  '/^  niul  the  homulary  conditions  lor  this  propciiv 
should  be  atlpisted  to  :ai^e  the  correct  are;i  under  llu'  curve. 

!n  tiddition  lo  featuring  the  ini  ra-elimient  modelinj’  ea[)abi  lily,  this  appli- 
cation also  rtiises  th«.-  issue  of  seiisitivilv  anahscs.  ,S<’hmil  h;is  Jong  advocated 
making  tmalysis  ]jr()grams  more*  design  t)riente<i  by  provitling  inloi  mation  on  the* 
sensiti\ii\  of  the*  imah'.-^is  rc,-.ults  to  (K.*sign  chtuiges.  I'lie  present  lormuhtlion  (d 
the*  noidinetar  male*ri;il  iltccls,  K((ualion  iMo.  allows  such  inlormation  Lo  ije  com- 
puU*d  by  making  the  change  in  m;ite  i*i;il  i)ropi  rlies  constant  with  re*s[)ect  to  a slate 
vtiriable  and  ludting  the  iteration  after  the  lirst  or  se'cmid  e\r!c.  This  procedure 
was  siK*cessfull\'  applied  to  a unit  cell  anal'.sis  as  pa  i*t  ol  another  studi'. 

1. 2 rrae*ke*d  liar 

d’o  establish  the  utility  of  the*  paramel  I'i/ation  in  Table  2 lor  tliree- 
dimensional  enick  anahsis,  an  internally  cracked  bar,  Figtirc  h,  studied  by 
Gyt'kenyesi  and  Mendel  son' “ was  analw'.ed  using  a two  element  semmelry  mcKlel. 
This  problem  is  a finite  dimensional  version  id  the  pentn  slui|)e*d  or  poker  chip 
crack  {problem  wliose*  elasticity  sobation  i-onUiins  a squtire  root  strain  singularity. 
The  two  element  PATCIIKS-Ill  model  shown  in  Figure  JO  uses  parametri/ations  in 
the  axial  and  radial  directions  that  induce  a scjuare  root  stnun  singularity.  Wlien 
the  singularity  is  at  Z (1)  nither  than  Z(o),  as  in  the  element  over  the  crack  face, 
the  formal  expression  for  /.  (fi  is 

7(0  a (21i 

which  lend  to  Ihe  s:ime  behavior  at  Z (M 

u,  j.  " 

-(a'2^'r)u,^  t‘-2) 

where  r l-7.j  in  this  element  and  the  comma  notation  indicates  differentation. 
There  ari*  two  fealiir(*s  of  Mu*  parametric  cubic  nuKlel  that  merit  attention  1.1 
the  gitometty  model  is  trivial  to  construct  in  algebraic  or  geometric  format  and 
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vai'i;il)le  whirl)  is  iliivi-lly  l•<-“luted  lo  the  stress  intensity  lac-tor.  Tl)e  eoeffieients 
for  the  al};el)i'aic'  , S,  , and  jteometrie,  Bj,  representations  of  L'quation  (21)  are 
si  inply 


and 


S 

1 


<h  s. 


-a,  S., 


, S , 


0 


(2d) 


Bj  U,  B^,  a,  B.j  2a,  B.,  - 0 (2-B 

The  eonstruetion  of  the  complete  jfeometry  model  for  this  problem  required  six 
)trid  cards,  two  PATCllB  directives  containing’  only  four  nonzero  B_  each,  and 
'wo  IIPH  directives,  liowevei',  this  easily  c-onstrueted  model  is  capable  of  accur- 
ately representing  the  iiighlt’  deformed  geometry,  Figure  11,  as  the  crack  opens 
under  an  axial  load. 


Comparison  between  the  present  solution  and  the  Gyekenyesi-Mendelson 

solution.  Figure  12,  show  good  stress  agreement.  The  strains  in  P.\TC  HKS-Ill 

are  computed  at  the  (;aussian  [loints  and  then  transformed  to  the  one-tliird  points 

which  accounts  for  the  finite  am|)litude  stresses  at  r = 0.  The  stress  intensity 

factor  computed  from  u,^(U  of  the  element  containing  the  crack  face  is  2.5'7 

higher  th;ui  th;it  for  an  inl'inito  dimension  bar  whch  is  slightly  closer  than  the 

(lyekenyesi-Mendelson  result.  The  displacements  are  also  correspondingly  closer 

to  classical  than  theirs.  The  parametric  derivative  u,^  in  the  present  model  may 

provide  a convenient  characterization  of  stress  intensity.  The  issue  would  depend 

on  its  behavior  under  admissible  repaitimelrizations  such  as  those  in  Table  2.  Since 

in  general  the  strength  of  the  singularity  is  unknown  and  changes  as  a function  of 

:?:i 

material  anisotropy  . A carjjet  plot  of  the  radial  strain  component  over  the 
(/^j.  /.j)  plane.  Figure  12,  shows  a theta  dependence  in  keeping  with  the  elasticity 
solution  which  contains  sin  0 and  cos  9 terms  in  the  strains. 


•I . C.raphite  Bar 

Consider  next  the  applit-ation  ot  a stress-strain  state  variable  model  to  a 
composite  material  whose  inelastic  behavior  is  known  to  differ  significantly  from  that 


i. 


- 21  - 


of  metals.  ATI-S  ^I'aphile  is  one  such  maliTial  and  its  inelasLie  bidiavior  has  heen 

•> 

extensively  tested  l>y  .loriner'’  in  biaxial  and  triaxial  stress  slates.  To  focus  i learlv 

on  the  abilitv  of  the  model  to  reptX'senl  inelasiic  behavior  not  normally  found  in 

metals,  a c:ise  of  hydiMstalic  compression  was  selected.  ,\  classical  [ilastieity 

solution  would  lie  identical  to  the  clasticilv  solution  for  hvdrostatic  loadinj;.  This 

material,  howiwer,  is  ciuiie  nonlinear  under  hvdrostatic  compression  and  tfie  Haldorf 

mode!  for  .VI'd-S  fits  the  -iornter  test  data  extremely  well,  l igaire  M.  A trans- 

versels  isotropic  };raphite  bar,  diaineti.‘r  ti.^aP  inches  by  l.u  inches  was  loaded  by 

an  axial  force  and  I'Xiernal  |)ressure  in  these  ti'sts.  The  stress-strain  version  ol 

the  Hatdort  mode!  v.tis  used  to  cominile  the  data  in  figure  1-!  assuming  constant 

strain  ratios  which  is  consistent  with  the  -loilner  results.  The  bar  elastic  constants 

and  I'  foi-  the  Batdorl  model  are  given  in  fable  7 where  the  across  grain  direction 
ijkl 

o o o I o ■ 

coincides  with  the  centerline  ol  the  bar  and  the  parameters  are  and  in  , 


These  cmisuints  are  determined  from  umcxial  d.ata  incliuting  the  constant  multiplier, 
f*’,  where  O.lb  for  .A  Ti-S.  J'he  use  of  constant  strain  ratios  in  the  a>;ial 
tension  range  causes  tlu‘  stress  ratios  t rom  the  Hatdort  moticl  to  deviate  sligiitly 
from  1 1 hut  this  difference  is  small  even  up  to  0.00-1  as  Figure  1-1  demonstrates, 
fable  7 (ira|)hite  bar  matt'rial  lonstants 
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A I'A'IVIIKS-lIl  Muxlul  ul  the  g'.ige  aeetion  of  the  test  speeimeii  was 
developed  iisin^  the  material  eonstaiils  Irom  Tiiljle  7.  I'he  {teometrie  simplieiiy 
and  transverse  isotmpv  allowed  a simple  one  element  symmetry  model  ol  a do*'’ 
sen'inenl  of  the  bar  to  be  used.  A sinuk'  hydiostatie  load  eondition  of  -1  ksi  com- 
pression was  analyzed  iisinj;'  successive  elastic  solutions  ol  liquation  (la)  starting 
from  the  linear  solution.  fonv(!rgenee  of  the  pseudo-loree  method  was  rapid  as 
Figure  la  illustrates  and  the  resulting  strains  are  in  good  agreement.  Figures  10 
and  17.  with  Jortner's  test  data.  In  tliis  particular  pi'oblem  the  exact  stress  solu- 
tion is  known  from  eqiulibrium  considerations  and  the  iteration  was  tei’minated  when 
the  stresses  weie  within  one  percent  of  the  exact  solution.  The  convergence  char- 
acteristics of  several  interesting  parameters  are  shown  in  Table  s. 

Table  s .Xonlmear  solution  convergence  data 


Cyle 

C 

Potential 

Knergy 

^2 

Q 

max 

0 

-4000 

-. 00014 

. 7537 

0.0 

1 

-374  8 

-.  00973 

. 5189 

2.33 

2 

-3845 

-.  00119 

. 0335 

3.  52 

3 

-3909 

-.00132 

. 0882 

4.  1 7 

4 

-3948 

-.  00139 

.7171 

4.  54 

5 

-3970 

-.  00143 

. 7331 

4.75 

GO 

-4000 

. 7337 

The  computational  performance  of  the  model  on  an  UNI  VAC  1108  was 
good,  taking  less  than  one  minute  per  pseudo-load  cycle  for  this  small  problem  and 
should  be  better  on  a parallel  processing  machine.  The  number  of  conjugate  gradient 
cycles  required  to  maintain  constant  solution  accuracy  decreased  with  each  pseudo- 
load cycle  but  not  dramatically.  The  original  linear  solution  starting  from  zero  used 
74  cycles  and  the  last  solution  required  54  cycles.  The  time  required  to  generate 
each  cycle  was  negligible  in  comparison  U)  the  matrix  solution  costs. 
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TliL'  system  desi^iieii  Uiid  dt‘V(.‘lo[)ed  in  this  stiidv  h.ts  ttie 

tle\il)ilitv  tt>  eneomptis.s  most  ol  the  available  inelaslie  eompusite  material  models 
and  appears  to  ln\  e the  ellieienev  to  make  its  use  leasible  il  nut  pitietieal.  It 
will  require  aiiditional  expeiieiiee  and  testing;  in  a vanetv  oJ  applieatiuns  before 
the  effieiency  of  the  system  ean  be  luily  established  but  the  basic  desi^pi  philo- 
sophy seems  sound;  namelt’,  better  mimputers  art-’  a more  likeh  de\elopment 
tlian  bette'r  numerical  methotis.  The  use  ol  stix'ss  or  strain  state  \aiiaJjle  modeiiiif^ 
is  certainly  not  new  but  the  de‘si;i,n  ol  one'  siaiulard  input  lormat  loi  a entiie  class 
of  malei’ial  models  is  a step  lor\sard. 

The  intiaxluetion  of  isoparametric  modeling  for  material  (iroperties  as 
well  as  geometry  is  another  important  step  required  by  the  variability  of  composite 
materials.  In  each  instance  mentioned,  the  results  presented  are  a beginning  with 
much  additional  uork  rtxjuircd  before  the  triaxial  inelastic  belia\ior  ol  composites 
is  encompassed  by  the  a\'ailable  models.  This  work  should  include; 

1.  Multiaxial  testing  programs,  such  as  dortners,  lor  other  three- 
dimensional  I'omposites  that  iiududes  a parallel  computational 
model  development  effort. 

2.  Kxtensions  to  include  interstitial  slip  with  friction  at  bimaterial 
interfaces.  'I'his  behavior  in  fiber  reinforced  composites  is  im- 
portant to  an  understanding  of  microst ructural  effects  on  the 
failure  of  these  materials. 

3.  IXivelopmcnt  of  alternatives  to  point  stresses  and  strains  as  a 
measure  of  composite  material  response  for  structural  appli- 
cations. 

1.  Development  of  representative  volume  element  models  suitable  for 
inelastic  behavior  and  failure  that  account  for  the  statistical  nature 
of  three-dimensional  composites. 
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i^a ranietric  Cubic  Models 
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MAXIMUAvA  DISPLACEMENT  IN 
SCALED  COORDINATES 


FiK^Te  0.  Energy  Convergence  of  the  Conjugate  Gradient  Solution 


Figure  7.  C a rlion- carbon  T'nit  cell  Schematic 


JORTNER  TRIAXIAL  DATA  

PATCHES  -m/BATDORF  El— EMU 


-0.001  -0.002  -0.003  -0.004 

ACROSS-GRAIN  STRAIN 


riRurt'  Hi. 


liielaslic  Across-grain  Strain  Convergence 


TRIAXIAL  STRESS  (KSI  ) 


-5 


Figure  17.  Inelastic  With-graiu  Strain  Convergence 
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Tht'  prugiam  updates  lor  nonlinear  maierial  nuMlelmti  in  l>A'I'C II KS-III 
were  made  usin^'  the  nenerali/.ed  poslpi'oeessor  system  (d‘()Sri>.  This  system 
interfaces  with  PATCllKS-lll  ihrousrh  three  files;  I’l'DATA,  INI>T  and  HANlXtMKi 
created  during  a normal  execution  of  a lineai'  elastic  analysis.  .Any  routine  in  the 
lAVrClIlhS-Ul  library  is  available  to  CIH  )ST1>  and  allows  new  capability,  like  .MA'l'N, 
to  be  I'ully  developed  and  tested  before  I'est ructuring  the  oiiginal  code. 

In  the  present  effort  a major  change  had  to  Ijc  made  to  the  code  to  cope 
with  large  nonlinear  matrix  problems.  Tliese  changes  are  complete  and  the  re- 
structured code  is  shown  in  the  update  t(j  Figure  1-2  of  the  programmers  manual. 
The  GP(.\STI>  system  for  the  M.VFX  postprocessor  is  shown  in  I’igure  l-2a.  The 
limited  testing  of  MATN  accomplished  dining  the  study  tuis  been  very  successful. 
This  -link  will  tic  adiied  to  the  basic  system  after  all  options  ai  e full}-  tested. 
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lnj)ut  Data  (Aird:  I'ATC'llGH  Patch  goiiorated  by  general  line  rotation. 


Deseription:  Generates  a bicubic  patch  for  the  surface  created  by  rotating 

a PC  line  about  a general  axis  of  rotation  through  gtinima 
degrees. 


Format  and  Example: 


PATCHGH 

ID 

n 

UD,  SKG 

ZAl 

PATCHGH 

5 

3 

1.5 

TID 

aAMMA 

CAM.\L\0 

■^Pl 

25 

l.ID,  SEG 


Contents 

The  identification  number  to  be  given  the  patch  generated  from 
line  LID,  segment  numljer  SEG. 

The  line  number,  l.ID,  and  segment  number,  SEG,  that  identifies  the 
PC  line  to  be  rotated.  A blank  SEG  defaults  to  one. 


/.'\[,  '/.BI  Coordinates  of  two  points  that  define  the  rotation  axis  directed  from 

2a  • 

I'lD  Transformation  ID,  if  any,  that  defines  a geometric  transformation 

to  be  applied  to  the  PC  line  before  rotation.  The  line,  LID,  does  not 
change. 

GAMMA,  The  angle  in  degrees  through  which  the  PC  line  is  rotated  starting  Yq 
degrees  from  the  initial  position  of  the  line.  The  sense  of  rotation  is 
determined  by  the  right-hand  rule  and  the  directed  line  (vector)  from 
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li^jut  n:ita  (\ii-d: 


I’A  I'CIK)  Oiilluiu  patcli(os) 


rv  s c ription: 


The  paleh(es)  >;eneralo(l  by  moving  an  ouUine  cun,'e  along  a 
base  eiirve  with  a fixed  orientation  of  the  outline  curve  in 
the  global  fniine  or  in  the  local  Frenet  frame  of  the  base. 


Forma  tt^'d  Jlixa  mplo: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

I>ATCtiU 

IDl 

Bun.  SKC. 

01.11),  pro 

TID 

y H:\yi' 

ID2 

IDS 

PATC  HO 

5 

6.1 

3 

F 

11 

1 II, iC 

15 

Field  Contents 

lUl  Patch  identification  number  for  first  patch  generated. 

BUD,  SEG  Baseline  identification  number,  and  segment  number,  SEG. 

If  the  SEG  is  not  specified,  the  entire  line  BLID  is  used. 


OLID,  SEG  Outline  cui-ve  identification  number,  and  segment  number,  SEG. 
If  the  SEC  is  not  specified,  the  entire  line  OLID  is  used. 


TID 


Transformation  ID,  if  any,  that  defines  a rotation  matrix  to 
reorient  the  outline  curve  relative  to  the  base  curve.  The 
outline  curve  is  always  translated  to  the  first  grid  point  of  the 
baseline  independent  of  TID. 


FRAME  E for  fixed  outline  orientation  with  respect  to  the  Cartesian  frame. 

p’  for  fixed  outline  orientation  with  respect  to  the  local  Frenet 
frame  of  the  base  curve. 


1D2,3, , N List  of  identification  numbers  to  bo  given  the  second  and  subsequent 
patches  generated,  if  any.  This  sequence  proceeds  from  the  second 
line  segment  of  OLID  to  last  and  then  repeats  from  the  first  segment 
for  the  next  segment  of  BLID. 
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miLK  DATA  I)i:CK 


Input  Curd:  MA  I'OH  Orlhotropic  Material  Definition 

IX’seription;  Defines  the  material  properties  for  a linear,  temperature 

independent,  orthotropie  material  from  engineering  constants 


Format  and  Fxa mple: 


2 

3 

4 

5 

6 

7 

8 

9 

10 

MATOR 

MID 

KHAME 

POINTS 

01 

02 

ON 

r>UTOR 

3 

1 

8 

2 

6 

THRU 

12 

Field 

MID 

FllAME 


POINTS 


Contents 

Material  identification  number 

= 1,  Pi'operties  are  in  an  orthonormal  Cartesian  frame 
-2,  Properties  arc  in  the  normalized  parametric  frame  of  the 
element.  Assumes  parametric  frame  is  quasi  cylindrical 
or  speherical. 

^1,  Constant  material  properties  (N=l). 

8,  Trilincar  variation  of  material  properties  (N=8). 

-04,  Tricubic  variation  of  material  properties  (N  d). 


01 ,02,  . . ,ON  Matrix  identification  number  for  the  orthotropic  material  engineering 
constants  at  the  inteipolation  points.  If  POINTS  is  equal  04  a 
single  entry  is  used  to  identify  a matrix  containing  the  04  OlD's. 


Roma  rks; 


1.  The  engineering  constants  are  entered  in  sequence 
^22-  E33,  Uj2,  0-^3,  U23,  G^2>  ^13,  C,23  on  the  MTHX 
matrix  card  (s). 

2.  The  Air  Force  Design  Guide  convention  for  Poisson  ratios 
is  used  i.e.  Ejj  i>ji  = Ejj  Ujj . 
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lUU.K  DATA  DIX'K 

Input  Data  Card:  MATN  Nonlinear  Material  Properties. 

Descriittion:  Detines  the  nonlinear  beliavior  ol  ;i  mtiterial  property  in 

ti'i’ins  of  stress  or  strain  sltile  varialjles. 


Format  and  Kxample: 


2 ;i  I 'i  '■  7 “ it  t'l 


M.\TN 

Mil) 

Ml’II) 

1 

M.AT 

SM.MM-; 

S'l.A'l!  1 

ST  ATI '2 



-J 

1 

■‘.Ml 

.MA  I N 

U) 

M.\T( ) 

IT 

. 

J 

•Ml 

■ M 1 

I'l 

.MTKX-M 

VI 

M l liX-Vl 

\-x 

'.n  KX-V2 

-M2 

1 

■tn 

dUO, 

ti 

IDO.U 

- M:j 

Fit-Id 


Contents 


MID 


Material  idcntilieation  number. 


MPID 


MAT 

SHAPE 

STATE I 
VI 

TI 

MATHX-M 

MTHX-VI 


Identification  number  of  the  matrix  that  defines  the  mesh 
points  which  have  these  properties.  The  default  is  all  mesh 
points  in  an  element. 

Mnentonic  that  defines  the  format  of  the  data  in  MTKX-M. 

Mnemonic  that  defines  the  interpolation  methotl  for  T,  \1,  V2 
in  that  order;  CEL,  l.CC,  etc.  where  L linear  and  C = Cubic. 

Mnemonic  that  identifies  state  variable  I (see  remarks). 

Value  of  the  variable  I for  which  these  data  apply.  A blank 
field  requires  .MTUX-VT  data  to  define  VI. 

Tempe-rature  at  which  these  data  apply.  Tlie  default  is  all 
temperatures. 

M.atrix  identification  number  of  properties  data  in  the 
format  defined  by  MAT. 

Matri.x  identification  number  of  coefficients  that  define  VI. 


Remarks 


1.  The  standard  state  variables  are  VI  and  \’BI  for 

I - l,2,;i,  l where  VI  = FjCj,  VBl  ^ F^jy  ‘^i°j'k'^l’ 

2.  llie  only  non-standard  option.s  are  the  Batdorf  model 
identified  by  a 13  and  the  Lee  identified  by  an  L. 
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Input  IXUu  Card:  TMOVK  Uigui  Body  Tran.st'ormalion 


Description:  Defines  a tratisfonmuion  that  move  objects  (lines,  patches,  inper- 

patches)  as  rigid  bodies. 


Format  and  Example: 

12  34567  S9  10 


1 r 

TMovt: 

IP 

1 

ZOI 

1 

Z02  1 Z03 

1 

DCIU 

V 

1 

e j 0 

1 

*TA  i 

m 

■■ 

3.  0 

0.  ! j, 

i 

30. 

20.  ' -10. 

1 

*T\ 

Tl 

T2 

1 

T3 

1 

! 

I 

3.  0 

0. 



4.D  ! 

1 

1 i 

i 

i 

Field 


ID 

ZOI 

DCID 


^,9,0 


TI 


Contents 

Transformation  identification  number  (1  to  100). 

Defines  an  origin  for  rotation  of  the  object. 

Direction  cosine  matri.x  identification  number.  If  blank  or  zero, 
the  Euler  angles  'f,  6,  p,  define  the  rotation  matrix. 

Euler  angles  in  the  3,  1,  3 rotation  sequence. 

(coi  ^ coi  ^ — CO#  9 ♦ sin  — sin  ^ cr»s^  - ro«  9 #in  * cos  ^ sin  9 #in^\ 

i-o#  ^ sin  I*  T ow  9 c«i#  « #m  ^ — sin  ^ sin*  ♦ t o#9  ccbj  * . i«s  — sin  9 cus*  J 

siD9sini^  sin9cosi^  co#9  / 

Defines  a translation  to  be  applied  after  the  rotation. 


He  marks; 


1.  The  complete  transformation  can  be  defined  as 


Carljon-C'ailjiMi  I'mt  C'i‘11 
Hulk  Data  M(xiul 


cr-rr-r-r-f'r^r'a  cc  a «: 

• Aj  A»  Aj 

oattcttctr.  era  xai/'CLraLPOir 


^ c c r c c rr  t 

rrcoe’-CF'cc 


. ru  — ^ A.'  •—  IT 


OOCOOOeOC  / 

« • > •••  •-A. 

c>  cr  c o c c-  c f. 


occ«~ooc  r 


OOOOOC.OC  C,  ' 

fi  K-  r^-  K'  f'-  f*  c- 
Cjccccccc;  • 

c.  C-  o c c.  c c c* 

f-  r-  r- 

iAipiTj  \r  s ^ 

ec-xaaaofT'i 

oocac  O O O r: 

Oa3nj3  J:J  *1 

o -c  -o  -c  C i c £ c:  > 

ooccotcc 

)COOCC  C'OC-OC 
,rsjrv,rv»r\jFAfo»^K>K'K% 


ococ-coor.  0( 


occ  c c oooc 


■ C l'*  c IT  c>  < 
o o c.  o o < 


c 0*000000000  cc;c-oo  i 
«>  «•  — ^ 

cj  c.'  o ^ ^ ^ -••  A.  3 ir  c-  I 

c CCC  C C C C C C-CC^CC-CCCC*-  >»,>«>*>».>«>.>' 

CD  ir  »/,  (TV.  viovr.  </»!/</  ir«/w?i/“i  ^ « 4^.^  4 ♦♦ 


».  r\.K  JJ-c^'C<^c■•-'A»'  J»^•C►-erco  — A;f^;3.r'C^  Tcrc— 

» f\,rvrv.fVAjfVA.r\.r\j«\,r»  r^'r--.  r*>r.rif 


PAG? 


^•<#7  7 


«-:fL'1-!,‘^Tf.l'>-?,t!FI.''-!,rTrur!-«,PTfL0-‘;,f!El.D-fr,FIEI,0-7,F!r'.  0-e,FIEl.D-'5,FTttD*10 


r 


— 


Oj 


>C  — 


^ a 


c:  — iT  -C 


^ — a 


IT  C- 

f\.  tr  ro 

cr  — - 


<£ 

C 

X 

r-- 

XI 

X> 

4^ 

r~ 

r- 

■C 

X 

■O 

X 

•C 

r- 

r- 

O 

o 

O 

o 

o 

c 

o 

C' 

o 

o 

o 

o 

o 

o 

o 

c 

c 

c> 

X 

o 

o 

c. 

o 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

•4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

;3 

iT 

K. 

f- 

IT 

• • 

d 

Tn. 

kT 

cy 

u 

C 

X 

X 

X 

X 

c. 

4’' 

o 

O 

-C 

r- 

IT 

C/ 

d 

a 

-■- 

IT 

U* 

C. 

cr 

«-<• 

d 

iT 

iT 

y 

If 

y 

<v 

d 

c 

(T- 

c 

r- 

<\- 

£ 

X 

o 

C 

c 

<r. 

■C 

<£ 

X 

a. 

cr 

a 

o 

cr 

d 

•— 

X 

a 

o 

a 

a. 

o 

— 

a 

4'* 

c.. 

4^ 

X 

X 

a 

X 

a 

— 

t~ 

■c 

r- 

r- 

r- 

0 

r 

K 

4- 

r- 

r 

c 

X 

X 

X 

X 

r- 

o 

o 

o 

o 

o 

c- 

r? 

o 

o 

C.' 

c- 

O 

kT‘ 

C' 

t- 

c 

o 

o 

o 

o 

o 

o 

o 

♦ 

♦ 

♦ 

♦ 

> 

♦ 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

c- 

o 

r- 

--• 

d 

C 

d 

o 

,c 

rv» 

a. 

X 

cr 

a' 

O 

y 

d 

a 

•1 

rvf 

lT» 

a 

d 

o 

»- 

IT> 

-C 

c 

r~ 

\r 

d 

ir 

tr 

IT 

•y 

y- 

f~- 

X 

c. 

w-t 

d 

X 

c 

r* 

a 

C 

f'l 

rv 

r- 

d 

p 

c 

cr 

o 

d 

X 

c 

c 

• 

— 

or 

- 

— 

o 

- 

'• 

•* 

-* 

X 

X 

a. 

ff 

X 

a 

c 

a 

C 

•c 

a 

c 

•t 

X 

X 

<X 

rx 

c 

r- 

X:  f- 

4 

■c 

fN.  C 

r 

C 

■C 

c 

r 

X 

C 

X 

c 

X 

r 

a 

c 

X 

C 

X 

X 

X 

r 

a 

X 

o 

c 

O 

o 

o 

.2 

o 

c.- 

C 

O' 

o 

C- 

o 

o 

C 

O 

c o 

o 

c- 

O' 

o 

c 

o 

c- 

<r 

C- 

<T‘ 

O' 

c 

c 

X 

r> 

o 

C‘ 

c. 

c 

V- 

c 

o 

♦ 

« 

♦ 

♦ 

♦ 

« 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

4 

♦ 

4 

4 

4 4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

c. 

c. 

•c 

X 

u. 

c- 

f- 

o- 

r- 

c 

r- 

u- 

rv 

'V 

X) 

X 

X ^ 

a 

rr, 

rv. 

»r 

— 

r»~ 

r\j 

XT' 

.J 

ir 

d 

IT 

d 

y 

-7 

y 

X 

X 

r\. 

X 

c_ 

d 

r- 

c 

c 

c 

c 

er 

r- 

r- 

X 

a 

r' 

3 

r~ 

X 

f'- 

rv 

ir 

c. 

4».  C 

O 

J" 

a 

rr- 

C 

f» 

X 

r»  t 

X 

4»> 

X 

r»* 

X 

o 

O 

c 

c 

C- 

r\. 

tr 

a 

o 

s 

IT 

a. 

o 

cC 

c- 

u. 

c. 

<L 

c 

cc 

A. 

-c 

fO 

— • 

d tx. 

c 

a. 

tr 

<v 

fT', 

d 

r- 

r- 

t-~ 

r- 

o 

y 

•3 

o 

y 

a 

y 

r\j 

rv 

rv. 

il 

rv 

IT 

- 

IT 

l/> 

w 

d 

ix 

4~- 

4- 

a 

d 

C 

— 

3 

" 

cr 

** 

rr-. 

K 1 

K' 

4^. 

r\. 

w 

f\/ 

d 

rv 

r\ 

y 

i* 

J"' 

<• 

r«- 

X 

cr 

c 

r\i 

4»« 

d 

w 

X“ 

r- 

X' 

o 

C 

f\. 

rr 

- 

y 

X 

o 

O 

K1 

c 

o 

O 

r^, 

Ci 

rj 

o 

d 

o 

o 

o 

J 

c 

J CT 

j 

o 

d 

c 

d 

\Z 

T 

C- 

3 

;• 

»“ 

iT 

c 

y* 

c* 

y 

c 

i/' 

c 

y 

c 

y 

c 

»- 

c 

c 

o 

c 

C- 

«-• 

c 

V 

c 

o 

» . o 

c. 

o 

— 

c. 

• 

X 

-— 

c 

»• 

C 

— • 

c 

X 

c 

• 

f* 

•- 

c 

•• 

c 

1 

1 

lT 

1 

•c 

» 

f“ 

t 

<r 

1 

c 

t 

o 

• 

•-. 

1 

F\ 

1 

1 d 

1 

IT 

1 

c 

1 

4~ 

• 

<t 

1 

c- 

1 

o 

1 

I 

rv 

1 

4»^ 

1 

1 

« 

y 

1 

X 

V 

f^. 

r^. 

X 

K" 

r»'. 

w 

r-- 

>• 

>■ 

d 

>» 

d 

y 

d 

M d 

■ta 

-■» 

>« 

d 

>* 

d 

>• 

^3 

m 

d 

>• 

y 

> 

J"' 

> 

y 

X 

y 

X 

y 

M 

y 

y 

rr 

•- 

a- 

or 

_< 

a 

fi- 

ur 

». 

a 

» 

ZI 

or 

a -• 

d 

d 

w- 

d 

a 

-•■• 

a 

•-• 

a 

a 

. 

a 

«• 

a 

»•> 

a 

•> 

a 

- 

.i 

-- 

» 

? 

p- 

3 

2 

► 

: 

»- 

• 

V 

► 

> 

4- 

1 

»> 

*- 

4 Z 

4 

3 

4- 

4 

3 

4 

J 

4> 

; 

4“ 

3 

4 

i 

4 

3 

4 

j 

4 

3 

4 

r 

4 

3 

1 

♦ 

a 

♦ 

♦ 

> 

« 

3 

♦ 

1 

♦ 

; 

* 

3 

4 

3 

4 

♦ 

5 4 

3 

4 

- 

4 

- 

4 

i 

4 

3 

4 

a 

♦ 

> 

4 

S 

4 

3 

4 

1 

4 

3 

4 

3 

4 

ac-»\»^^7'£^c/<7c- 


• 'V'orjtT  7 

4/^,Ti/  irC-CCA4  •>.  c c *.  L 


- fv  f->  ) O’-  c 

f'  r-  ^ r f- 


a a ^ 

r-  r-  .i  .1  •»  7 «*. 


\ 


CAR?  crc,__n-!,rtfLr-?,F!FLn-!,CTr|o-u,rTtL0»S,FlELn-6,FIEL''-7,FTrLn«6»FIEL';->’,F!F.LO»10 


ir'iririr^u'*rw'»rirj^ir'ir‘iri/'*rkrj"j^'»riririr^»ri/"*rir.r»ir 
oc-cc-c'c.cc  c-cooocrc  C.OSC.C  ccooccroc-j'o 


— r^K".  o — ooa  f^—  r-.  i/r^r'ia.'^>afV»-a  ■cci/' 

r-K.  Ci;rv.<\,ia  oaaa  — iTcc''  — « 

iT-ocj"kr*(ycr  '\jyxj“xJ^it»-3  j'c.'v.'^r*.  c c oci 


X X r-  r-  -o  X r- 

X o c.  o c o x*  c crxrL'Cor  =cox'C*rix-xc  _x-r  x-c: 

♦ ♦ ■¥  ♦ ♦ ♦ * * •lfl|llllllll«llllllllllllllll 

X •-  X o c a X rv,  r-xacr*-r*‘r~r^r^tirc*r-i/**'— 

cj  r\.  a — c.  ^ X »»  :5t^eacr-*'c^—»*  -.*c,^r-i/-rvjr'crrcccu-  — XX 

iT  r-  o f\j  c-  X — ca'’^j'cc.x.  cx\.coxoa^-r^c--xr\.^ccx''vt-*T 

r-  (t,  — — r-  cc  •“  ^ - rv 


r-  xr~  xr^  xr~  X c x xx»^  x^-  xu^i/  J'lTir.xirj^iXj'j'Xirirj  j~  s o J"  J-  >r  iT 

cor>o'TO'ic>';^:o*roc:wC-:.cc-c-CiC'>crr".  ccc-'i^oxcc'X-  OsTcc  r t.  c> 

fv  -•  J a r . r?  C7  K^I  C X rvt  IT  — C C — u*"  rv  -'•  J — C — C"  X C »r  — — — 3 

X^o-J^crv^*  *7C*J*^C  xu'^cr*'^»'cf*-r-»'‘r*f^r'r-f^r^r'r*r^*“rv»a^c‘-i/'ccxcr^r' 

► rv.'wxx'r  *.ooc  xx-f^  ^-iiXcOi/j'—  — r*'0  c"^c  xccaac-  — 

O O jryXOf^fr^ffCJ  — r^  - — 


r <x  o c •“  <v  -■»  ^fV*^'ryxr'ao‘C^rv^-;3i/xr-i.ac-rv^rArxr't(7  '- 

irxt'‘oj'occ  xc’xcccrcooccococr:  — ^ — — • - — rvr^rvv  »vrv  r^.x'V'V'^ 

^o— C*  — o-C-»  c‘.-or\jrvxr\ir\.r\jr\,r\ir\r\.»\,'vrvrs.rvrvrvr\,'V^'V'V»\,'X'\»r\v<\rv,fv 

>«k/’>'l/'>«tX>*  X K X M 4 >»X*»  X»^>*>»>'XX>*>>»»»'>«>>*»*>«»*>*  XX  X ^XXXXXXX 

rr»-ri»  r*-3rr>T'i'iii  r^ir3'xr»af»ir'if»»-3r>a 

X Jx  ix  2x-3»  : ► :►  !►  ►►►  »■  K X ► x»  »■  x» 

1 *3  *3  *3  *^3*3  *3  *3.31313  33  3.33Z3  3X3  33  3.33  13  3 33  ^3  33 


iX  f X-  1 cr  o X <v  f ^ u (?  c<  fv.  X T j"-  X a C c.  — rv  X c r-  u.  o - rv  - .j  ^ ’ c 

, ..<V'\ifv'V'VrvrV  \ rvr\,xi 

x«»*«..  ..^•.x-«^x«*  •••x-fvrvrvjrvrvrvrvrvrvr‘»>/\jrvrvrvrvrvrvrvrvr\,f\.r^fvrv'V<V‘Vr^rwV<\j 


F’r',,')-’.  ,i!Fir-2<‘-'!FLP-3<ftfLC*‘<'.'='rf.L0-?iFrFL5-fc.fTE'l'5-7,Frt'. 


r 


■« 


o 


o- 

CJ 

u: 


CL 

CJ 


t 


o 


i\j 

t 

(V  ir  -c  o a.  IT  ^ *-* 

o 


O <£  r-  r\j 


^ — rv. 


•—  IT'  w—  •-• 


ir 

J' 

ir> 

IT 

J 

-.T 

iT 

J' 

IT 

tP 

ir 

IT 

\r 

uT 

4T 

vT 

ir. 

X- 

y 

IT 

tr 

o 

c 

o 

o 

C.' 

c 

c> 

c. 

O 

<_' 

»- 

c 

c: 

cr 

c- 

O 

' •' 

• 

1 

1 

1 

1 

ru 

A 

f- 

»>». 

r- 

r. 

<7 

c 

<7 

o 

r\j 

V/- 

(X. 

P- 

K . 

C7 

S' 

O <7 

X 

rvj  — • 

rv. 

C 

tP 

<v 

a 

ct 

r^\ 

r- 

c 

*3 

* 

CT 

•3 

X 

c. 

C7 

X 

X 

rL 

F\j 

-> 

f\. 

X 

C 

X 

c 

-C 

— • 

- • 

<7 

c; 

r\d 

rsj 

— 

- 

r~ 

rv. 

J 

*o 

*3 

.1 

o. 

Ts. 

t* 

rv. 

r\ 

* 

\r 

J' 

ir 

CJ 

c; 

o 

IT 

j* 

1/ 

\J 

J 

\T 

lT 

»/ 

S 

v/ 

J" 

J- 

L' 

J- 

\r 

J* 

v'’ 

X 

vT 

ir 

J' 

iT 

vT 

<r 

1/ 

<s 

LT 

y 

o 

o 

- 

-• 

- 

C 

O 

Cj 

o 

C 

' 

Y 

c 

- 

"-J 

'J 

Y 

J 

K" 

.J 

o 

»<■> 

r-- 

J 

c 

c 

J 

c 

- 

c 

c 

f' 

KV 

— 

C 

rL 

?7 

X 

— 

J 

pv 

.-3 

T 

r 

c 

o 

C7 

i 

X 

' ' 

*" 

■*" 

’■* 

^ 

f 

<7 

'L 

r- 

r\» 

J 

'V 

X 

C7 

!/■ 

- 

•— 

•— 

— 

:j 

X- 

’w 

t7 

— 

-J 

c; 

**• 

K. 

J 

a 

rv. 

U 

<\. 

'V 

iT 

J- 

vT 

O’ 

•3 

;3 

3 

iX 

' 

* 

ir 

•T 

J' 

\S 

J" 

J 

iT 

vT 

IT 

X 

X 

X 

■C. 

iP 

vT 

L'' 

y 

J- 

ir 

s 

v/~ 

vT  -r 

lT 

y 

o 

c 

o 

o 

•3 

c- 

c 

O 

o 

- 

~ 

- 

= 

"= 

o 

o 

y 

o 

C 

1 

xTl 

t' 

ro 

K* 

a. 

a 

X 

O' 

iT 

KV 

t7 

pr> 

ro 

PP  rv,  rv. 

-3 

ij 

;> 

J- 

to 

r « 

r- 

K' 

f*-- 

rv. 

X 

f\. 

X 

X 

a 

rv. 

— 

c. 

t" 

— 

o 

<7 

a 

3* 

<7 

»»- 

■c 

c> 

a 

»- 

%r 

•- 

— 

•— 

•— 

o 

T 

ir 

o 

V/ 

— 

=» 

- 

" 

- 

r\* 

r\_ 

K 

r<"i 

3 

J” 

J' 

sf- 

if 

• 

• 

• r\K*:j^^r'a  c ^ r\.  -c  ft  o c 

j - --J  r?  CT  r '•y' 

• I I t I f I t I I I t ( I I • * * * ' 

>fwn  yr  -»$  >•>*>'-"  M >« 


-rv^;JJ“  ^ ^ O C — ^ 

C € -T 

'\/'\jfV<\.r\.rvcv<rv.rLfL  ^'V'V  cv 

||*|*|llllll»»0000 

a rr  n a a jr  a :x  3-  ^ - 'i  ^ L.  LJ  c_ 

^ »-►-►  ►CC'*'  f' 

a.  3 3 2 3 3 3 3 3 1 3 3 a 3 U,  C'.  C'.  U. 


LJ 

lO 


- • a.  X 


rv,  rv.  r\j 


rL  o'  f 


t k/  s <j  c - 'v 


- f)(i  - 


y ■»■ 


* 


A 


s'^r  1 

SPC  1 0 

s r '.  0 


Od  :3  IT  — 


fV  n T IP 


t- 


« ;r  J ^ Of’  o 


^ :T  - — X> 


^ C IT  <£> 


• rv. 


rr\  1^,  f^,  r^t  r\t  e\j  >\.  fv.  Tv.  <X 


<^OCO  O'T-OOC'  c 

rrrir*'.  rrocrcc 
ir.  ir>  %r.  kT.  ^ ^ ‘'J 


^~f•^'tt.l.a  4.a  .\aa‘i  '-*- 


- ;>/  - 


^croe  - 


A^cpt 


l{Kl>tJKT  DISTUllil  riON 


FAirr  1 - c;o\  krnmknt 

Aili\uni^li  ativc  L.iiusoi^  Activities 

Chief  of  Nav;iJ  Itescareli 
IX'pnrtment  of  the  Navy 
Ai'lington,  Vir>;inia  22217 
AUii;  Code  -174  (2) 

171 

222 

Assistant  Chief  for  'reehnolo;2v 
Office  of  Naval  Research,  Code  200 
Arliiu;ton,  N’irj^inia  22217 

Director 

ONR  Branch  Office 
•1!)5  Summer  Street 
Boston,  Massachusetts  02210 

Di  recto  r 

ONR  Branch  (Office 
219  S.  Dearborn  Street 
Chicago,  Illinois  GOOO-1 

Director 

Naval  Research  L;iboratory 
Attn:  Code  2(i29  (ONRL) 
Washington,  D.  C.  20290  (G) 

U.  S.  Nav;il  Research  Laboratory 
Attn:  Code  2()27 
Washington,  D.  C.  20290 

Director 

ONR  - New  York  Area  Office 
712  BroacKvay  - 2th  Floor 
New  York,  N.  4'.  10002 

Di  rect(>r 

ONR  Branch  Office 
1020  K.  Oreen  Street 
Pas:iflena,  Calihjrnia  91101 


IX'fense  Documentation  Center 
Camei'on  Station 

Alexandria,  \ irgiiua  2231!  (12) 

Annv 

Commanding  Officer 

L'.S.  Army  Reseai'ch  Office  IXirham 

Attn:  Mr.  J.  J.  Murray 

CRD-AA-IP 

Box  CM,  Dul-ic  Station 

Dui’ham,  North  Carolina  2770G  (2) 

C o m m a nil  i ng  Off  i c -e  r 

AM.M\MR-ATL 

Attn:  -Mr.  R.  Shea 

I'.S.  Army  .Matei  ials  Res.  .•\gency 

Watertown,  Massachusetts  02172 

Watervliet  Arsen.al 
MAGGS  Research  Center 
Watervliet,  New  York  121 '9 
•\ttn:  Director  of  Research 

Technicid  Library 

Redstone  Scientific  Info.  Center 
Chief,  Document  Section 
U.S.  .Army  Missile  Command 
Redstone  Arsenal,  Alabama  35b09 

Army  R&D  Center 

Fort  Belvoir,  Virginia  220G0 

Navy 

Commanding  (Officer  and  Director 

Naval  Shi()  Research  & Development  Center 

Bethesda,  .Maryland  2002-1 

Attn:  Code  012  (Tech.  IJb.  Br.  ) 

17  (Struc.  .Mech.  Lab.  i 
172 
171 
177 

IsOO  (.-\(ipl.  Math.  laib.  ) 


-1- 


C'ommaiuiin^;  oHic  cr  ami  Dirt'cltji' 

Naval  S1h|)  lit-sfa la-i)  \ Ix-VL-Iopment  Cetitoi- 
Hethcsila,  .Marylaiul 

Attn:  C(Hle  n tliiS  (l)i-.  W,  1).  Settc) 
l!t  (1)1.  M.  M.  Sevik) 
l!H)l  (Dr.  M.  Slrassbcrg) 

1!)  ir, 

(Dr.  D.  Foil) 

1!)(I2 

Naval  Woapoiis  Laboratory 
Dahl>;rc'n,  Vii';j;inia  22Ms 

Naval  Hcseareh  Lalioratory 
\\'ashin<>  ton,  D.  C.  20;S75 
Attn:  Code  S-lOO 
s 110 

slljo 

SI40 

0300 

03;)0 

0330 

Fndersea  Kxplosion  He-eareh  Division 
Naval  Ship  K&D  Center 
Norfolk  Naval  Shipyard 
Flirts  mouth,  Vii'f>ini;i  23709 

Attn;  Dr.  K.  Ihdmer  - Code  7m0 

Naval  Shiji  Hcseareh  & IX-velo(imenl  Center 
/Knnapolis  Division 
i\nnapolis,  Maryland  21102 

Attn;  Code  27 10  - Dr.  Y.  F.  Wan^ 

2s  - M:  . H.  .1.  Wolfe 
2sl  - Ml-.  H.  If.  Niederherjier 
23M  - Dr.  11.  Vanderveldl 

Technical  Library 
Naval  Underwater  V\eapons  Center 
Fasadena  Annex 
3202  K.  FoothiU  Blvd. 

Fasadena,  ('alifornia  91107 

U.  S.  Naval  Wea(ions  Center 
China  l.;ike,  California  DSOr)? 

Attn:  Code  1002  - Mr,  W.  Wcrback 
1320  - Mr.  Ken  Bischel 


Comnuuuling  Officer 
( . S.  Naval  Civil  Fngr.  I^ib. 

Code  17)1 

Fort  llueneme,  California  1)3011 

Technical  Director 

U.  S.  Naval  ( )rdnance  laiboratory 

Wliite  Oak 

Silvci'  .spring,  .Maryland  20910 

Technical  Director 

Naval  Undersea  H<1;U  Center 

.San  Diego,  Ckdilornia  92132 

Supervisor  of  Shipbuilding 
U.S.  Navy 

Newport  New  s,  \'irginia  23607 

Technical  Director 
■Mare  Lsland  .Naval  Shipyai'd 
N’allejo,  California  94592 

U.S.  Navy  L’nderwater  Sound  Ref.  L.ab 
Office  of  Naval  Research 
F.  O.  Box  3337 
OiTando,  Florida  32306 

Chief  of  Naval  Operations 
Dept,  of  the  Navy 
V\  ashington,  D.  C.  20350 
.Attn;  Code  Op07T 

Strategic  Systems  Project  Office 
lX'[);u  tment  of  the  Navy 
Wa.shington,  D.  C.  20390 

.Attn-  NSF-001  Chief  Seientist 

Deep  Sul)mergence  System.s 
Navid  Ship  Systems  Command 
Code  39522 

lX)pm'tment  of  the  N;ivy 
Washington,  D.  C.  20360 

Kngineei'ing  Dept. 

U.  S.  Naval  Academy 
Annapolis,  Marykmd  21102 
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Naval  Air  Syslc-ms  Cummaml 
l)t'|)l.  of  Iho  Navy 
\\'ashinf;ton,  I).  t'.  2U.'iiiO 

Attn:  NA\'AII{  .».!U2  Aoro  ii:  Structures 
5;!(is  .sti'uetures 
r)2o;ilF  Materials 
CiO-l  Tech,  library 
;!2ol?  Structures 


.Ml’  l''orce 

C'oiiunander  W ADI) 

W j'if^ht-Pattei’soii  Aii'  Force  JJase 
Dayton,  Ohio  l.^Fid 
Attn:  ('(xle  WAVR.MDD 

AFFDI.  (FDDS) 
Structures  Division 


Director,  Aero  Mechanics 
Naval  Air  Develoiiment  Center 
Johnsville 

Warminster,  Pennsylvania  1897  1 

Naval  Facilities  Fiigineering  Command 
Dept,  of  tlie  Navy 
Washington,  D.C.  2U;J(iO 

.'\ttn:  N.WF.AC  03  Research  & Dev. 

01  Research  tv  Dev. 

1-4111  Tech.  Library- 

Naval  Sea  Systems  Command 
De[)t.  of  the  Navy 
Washington,  D.C.  203(i0 

Attn:  NAVSHIP  03  Res.  & Technologj' 
031  Ch.  Scientist  for  R& 
03-112  Hydromechanics 
037  Ship  Silencing  D'v. 
03.5  Weapons  Dymamics 

Naval  Ship  Engineering  Center 
Prince  George's  Pla/.a 
Hyattsville,  .Maryland  207s2 

.Attn:  NAVSEC  0100  Ship  Sys  Engr  & De:- 
0102C  Computer-Aided  SI 


AFLC  (.MCEEA) 

Chief,  /\{jplied  Mechanics  Group 
I'.S,  Air  Force  Inst,  of  Tech. 

Wright- Patterson  Air  Foi’ce  Base 
Dayton,  Ohio  1.5133 

Chief,  Civil  Engineering  Branch 
W l.RC,  Iteseai  ch  Division 
Air  Force  Weapons  Laboratory 
Kirtland  AFB,  .^e\^  .Me.xico  H7117 

.Air  Force  Office  of  Scientific  Research 
1-100  Wilson  Blvd. 

Arlington,  \'iiginia  22209 
.Attn:  Mechanics  Div. 

NASA 

Structures  Rese:u’ch  Division 
National  Aeronautics  & Space  Admin. 
Langley  Research  Center 
Langley  Station 
ILmipton,  Vii’ginia  23305 

IXjp 
ip  Des 


010.5G 

0110  Ship  Concept  Design 
(il20  Hull  Div. 

01201)  Hull  Div. 

012s  Soi  l ace  Sliip  .struct. 
0129  Submarine  Str  uct. 


National  Aeronautic  & Space  Admin. 
Associate  Administrator  for  Advanced 
Researih  K-  Technology 
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Nonlineai'  Composite  Materials  Solid  Geometry  Modelinf;' 

Triaxial  Sti’oss  Analysis 
3D  Isoparametric  Fitiite  Flements 
Carl)on-Cai’bon  Composites 

aBSTJUCT  'Conr/nut*  on  r«*verie  iide  /f  neres*ar>'  and  /denfifv-  hv  blocM  nufrbe'' 

A computational  model  for  the  analysis  of  structural  and  microstructural 
behavior  in  general  solids  of  composite  material  is  presented.  Fmphasis  is  placed 
on  reprcsentitig  tlu;  luiomalous  material  behavior  of  composites  and  on  the  construc- 
tion of  computational  nuxlels  with  vari;iljle  properties.  Alternative  maU>rial  models 
using  continuum  and  statistical  mechanics  were  reviewed  and  a modular  code 
designed  for  compatibility  w itli  several  different  models.  The  constituent  materials 


' ai  o l.•hal■actc'^i/od  in  terms  ol  those  st;ite  valuables  that  eorielate  a matei'ials 
response  sueh  as  ell'i'etive  strc'ss  oi'  stiaiin  enei’ity.  A parameti  ie  eubie  repi’e- 
si'iUation  is  used  loi'  all  state  variables,  the  solid  geometi'y  and  all  physieal 
(jroperlic's.  The  assoeiated  I'iidte  element  extends  iso|)ai  ametrie  modeling  to 
allou  piopei'ties,  linear  or  notdinear,  to  vary  over  the  volume  <d'  an  element  as 
in  rosette  material  eonstruidion.  Applications  o!  the  moilel  to  a earbon-earbcjn 
unit  ci'll,  t(,  strain  singularities  and  to  the  inelastic  response  of  a graphite  bai' 
illustrate  its  utility,  (lood  agrei-ment  \v  ith  triirxial  test  data  loi'  inelastic  sti  ains 
under  hydrostatic  pressure  is  obtained.  .Numerical  results  are  com|Hited  using 
I’A'l  t'l 1 1>-  Hi  and  the  conjugate  gradient  algorithm  without  the  assembly  of  large 
matrices.  This  approach  is  tailored  for  vector  processors  and  can  I'eduee  the 
high  cost  of  nonliiu'ar  three-dimensional  analyses. 


